Two Anderson impurities in a 2D host with Rashba spin-orbit interaction 
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We have studied the two-dimensional two-impurity Anderson model with additional Rashba spin- 
orbit interaction by means of the modified perturbation theory. The impurity Green's functions 
we have constructed exactly reproduce the first four spectral moments. We discuss the height and 
the width of the even/odd Kondo peaks as functions of the inter-impurity distance and the Rashba 
energy Er (the strength of the Rashba spin-orbit interaction). For small impurity separations the 
Kondo temperature shows a non-monotonic dependence on Er being different in the even and the 
odd channel. We predict that the Kondo temperature has only almost linear dependence on Er 
and not an exponential increase with Er. 

PACS numbers: 71.70.Ej, 72.10.Fk, 73.20.At 



I. INTRODUCTION 

Recent advancement in spintronics [l|, Q results from 
an increasing interest in studying the spin-polarized 
transport through nanostructures Q. Significant part 
of the research effort has been focused on the physical 
consequences of the spin-orbit interaction (SOI) and, in 
particular, on its effect in electron transport. SOI is a 
relativistic effect which is manifested in a coupling be- 
tween electron's orbital motion and its spin [J]. There 
are various proposals for spintronic devices in which the 
SOI plays in important role [2|. 

Kondo effect is one of the best-studied examples of 
correlation-driven phenomena in condensed matter sys- 
tems. It consists in the screening of the magnetic moment 
of a localized impurity by the Fermi sea of itinerant elec- 
trons [5] . Technological advancement allowed the Kondo 
effect to be studied experimentally in quantum dots J6j 
and in magnetic atoms placed on a metallic surface [7|. 
In these systems the delocalized electrons are confined to 
move in two dimensions. The inversion symmetry is bro- 
ken by the confining potential and the so-called Rashba 
spin-orbit interaction (RSOI) gives rise to a coupling be- 
tween electron's spin and momentum 0, Also, re- 
cently, investigation of magnetic impurities on graphene 
has been initiated [9(. Thus, a problem is posed: what 
is the influence of the RSOI on the Kondo effect in 2D 
system? Several studies have addressed it and the work 
has been centered on the Kondo model. It is obtained 
from the generic Anderson model by a suitable transfor- 
mation which integrates out the charge degrees of free- 
dom for sufficiently large Coulomb interaction between 
two electrons on the impurity |5|. In Ref. 10 an ef- 
fective two-channel Kondo model is derived and is pre- 
dicted that the RSOI does not change the Kondo tem- 
perature except for some band-width effects. In a later 
work [11], an effective two-channel Kondo model with 
additional Dzyaloshinsky-Moriya term is derived. It is 
argued that the latter interaction produces an exponen- 
tial increase of the Kondo temperature. Recently, Zitko 
and Bonca [Hj on the basis of an effective single-band 
impurity model concluded that RSOI leads to only small 



variation of the Kondo temperature with the strength of 
the RSOI. Thus, the behavior of the 2D single-impurity 
Anderson (or Kondo) model with RSOI is not yet fully 
understood. 

On the other hand, significantly less work has been 
done on the two-impurity Anderson model in a 2D host 
with RSOI. The two-impurity Anderson model (TIAM) is 
the simplest model that describes the interplay between 
the electron correlations and electron coherence which is 
the basic physics to be inferred from it. There are differ- 
ent regimes in the model due to the competition between 
the direct Kondo interaction and the indirect exchange 
(the so-called RKKY interaction). The former gives rise 
to the single-impurity Kondo effect while the latter is 
responsible for the relative orientation of the impurities' 
spins. Various approaches have been applied to study the 
TfAM - renormalization group approach 1 1311 . numerical 
computations T4J, perturbative method [la ], and modi- 
fied perturbation theory [l6| . Recently, study of the in- 
fluence of the spin-orbit interaction (of both Rashba and 
Dresselhaus type) on the direct Kondo and RKKY in- 
teractions has been performed in the framework of the 
two-impurity Kondo model [17| . 

In the present work we study the two-impurity An- 
derson model in a 2D host with RSOI. We use the so- 
called modified perturbation theory which has been suc- 
cessfully applied to both the single-impurity [l8j and the 
two-impurity Anderson model [161 ] . In this method an ap- 
proximate electron interaction self-energy is constructed 
which interpolates between the limits of strong and weak 
coupling to the band electrons coinciding with the exact 
expressions in both limits. More specifically, we con- 
struct electron Green's functions that are exact up to 



order 



This approach is not fully exact but gives 



very good results for only moderate numerical work. 

The paper is organized as follows. Section fl intro- 
duces the two-impurity Anderson model in a 2D system 
with Rashba spin-orbit interaction. The host part of the 
Hamiltonian is diagonalized and the corresponding quasi- 
particle energies are obtained. Section III is devoted to 
the construction of the impurity Green's functions in the 
framework of the modified perturbation theory and in 
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Section IV the results of the paper are summarized. 



II. THE MODEL HAMILTONIAN 

The Hamiltonian of two Anderson impurities, coupled 
to a 2D system with RSOI, is given by H = Htiam + 
Hrsoi with Htiam being the Hamiltonian of the two- 
impurity Anderson model 
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Here, Ci a is the annihilation operator for an electron 
with spin a =j, J, and energy residing on an impu- 
rity i — 1,2. ako- is the annihilation operator for a band 
electron with wave vector k and energy tk- Ui is the 
Coulomb repulsion energy between two electrons with 
opposite spins simultaneously occupying given impurity 
and V is the coupling between the impurities' and band 
electrons. The particle-number operator is = cf^c^. 
The impurities are located in positions Ri with the dis- 
tance between them being R = |Ri — R2I. Throughout 
the paper we consider only the case of two identical im- 
purities - 61 = €2 = e and Ui = U2 = U. 

The Rashba spin-orbit interaction is described by the 
Hamiltonian 

Hrsoi = <xr / J (fey + ik x )a£ f aki + H.c. 

k 

= a R J2(ke- i ' Ph a+ t a ki + H.c.) (2) 

k 

where cxr is its strength and the phase tpk is defined as 
k x = -ks\iiip k ,k y = cos tfk- 

In the presence of the RSOI one has to introduce the 
angular momentum basis for the band electrons. In po- 
lar representation the band electrons operators are de- 
veloped as 
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Ho is the Hamiltonian of the decoupled impurities' and 
band electrons 

H = n pa + — (rvr + n ot )(n ei + n oi ) 



pa 

E 

rah 



( m+l/2\ + m+1/2,, 

e k h [a kh ' ) a kh dk. 



(5) 



Here, p — e/o and ekh = £fc + cxRhk are the renormalizcd 
band energies. The term describing the hybridization 
between the impurities and the band electrons is cast 
into the form 
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where /3 e/om = (1 ±{-l) m )/y/2 and J m (z) is the m th 
order Bessel function |19j]. Let us point out that unlike the 
case of a single Anderson impurity the impurities' states 
couple to all orbital channels for the band electrons. In 
the former case, the impurity states couple only to the 
m = channel. 



III. THE MODIFIED PERTURBATION 
THEORY 



In order to discuss the dynamics of the impurities' 
electrons we introduce the retarded Green's functions 
G pa {t) = — i6(t) ({c ptr (0), c+ (i)}) (the curly brackets 
denote the anticommutator) . Note that the impurity 
Green's functions are diagonal in the even/odd basis. 
The Fourier transforms of the retarded Green's functions 
define the interaction self-energies E™*(w) for the impu- 
rities' electrons 

G pa ( W ) = [w - e - E<°? M - E™' ( W )] - 1 . (7) 
Ep 1 ^ (w) is the Hartree-Fock elastic self-energy [2(| 



where m is the orbital magnetic quantum number. The 
Hamiltonian for the band electrons is given by the sum of 
the third term in Eq. (JlJ and the RSOI. It is diagonalized 
by the canonical transformation 
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which introduces the chirality quantum number h = ±1. 
Next, we switch to the even/odd basis states for the 
impurity electrons (in the case of identical impurities) 
c e/ocr = (cict ± c 2a )/V2. With all these transformations 
the Hamiltonian of the system becomes H = Hq + H^ybr- 



E 



l±J (kR) kdk 

J - £kh + i0 + 4-7T 



(8) 



To proceed with further discussion of the physics of 
the system we are studying we must make some kind of 
approximation for the interaction self-energies. In this 
work, we choose to employ the so-called modified pertur- 
bation theory. In it the interaction self-energy is approx- 
imated by the ansatz 



E^*(o;) = -<n e _ CT + n _ a } + 



(9) 
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The first term is the contribution of the Hartree-Fock de- 
coupling of the Coulomb interaction term, Yy p J(uj) is the 
second-order self-energy which is obtained in a perturba- 
tion theory with respect to the Coulomb repulsion energy 
U, and the coefficients A a and B a are to be determined. 

(2) 

The explicit expression for Yr pa • (w) is as follows 
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The auxiliary quantities are defined as F(uji, u>2, ^3) = 
[1 - /(wi)]/(wa)/(u 3 ) + - /(w 2 )][l - /(W3)], 



5i(wi,w 2 ) = pi 0) (wi) / 9e U ''(w2) + /9o U,l (wi)/Oo Uj (w2), and 



5 , 2 (wi,w 2 ) = pi°''(wi)po U; (w2) + po U; (wi) / 9e U ' l (a;2). Here, 
/(w) is the Fermi-Dirac distribution function and 

Pp ^ (oj) = — Im Gpa (u) is the Hartree-Fock spectral func- 
tion. In the modified perturbation theory, the second- 
order self-energy is calculated using the following expres- 
sion for the Hartree-Fock Green's functions 

G p °J(uj) = [uo-~e-U{n e - a +n - a ) /2-E&>(a;)]-\ (11) 

the impurity energy level e being replaced by an auxil- 
iary energy level e which will take into account various 
renormalizations due to the Coulomb interaction and the 
charge transfer between the impurities and the bands. 
Later on we shall briefly discuss the ways to determine e. 

The electron Green's functions given by Eq. ([7} 
with the interacting self-energy from Eq. ^ may 
be constructed in such a way that the first four mo- 
ments of the corresponding spectral functions are exactly 
reproduced |2l|, that is, the G pa (uj) we are going to ob- 
tain will be exact up to order cj~ 4 . To that goal, we 
develop the exact Green's function G p<J (u>) [Eq. (J7J] in 
series with respect to 1/uj 
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The exact expression for the moment M pa ^ is given by 



MW = ({L% ff) c+}) (13) 

where the operator L acts as Lc pa = [c pcn H] (square 
brackets denote commutator) [2lj . Next, we develop 
the approximate Green's function G pt7 {uj) given by Eqs. 
(|7I9I10|) in series with respect to 1/uj. This is achieved 
by developing the explicit expressions for the elastic self- 
energy Sp5(w) [Eq. ([8])] and the second-order self-energy 
Yr p J{uj) [Eq. ([TO])] in series with respect to 1/uj. We are 
able to determine the coefficients A a and B a in Eq. © 
in such a way that the first four terms in Eq. (fT2")l co- 
incide with the first four terms in the expansion of the 



approximate Green's function, that is the approximate 
Green's function reproduces exactly the first four mo- 

(n) 

ments M pa ,n — 0,1,2,3. Let us point out that the 
second-order perturbation theory with respect to U (it 
corresponds to A a = 1 and B a = 0) reproduces only the 
first two moments. We shall not give details of the calcu- 
lations because they are too cumbersome to be presented 
here. The calculations along these lines but in the case 
of the infinite-dimensional Hubbard model can be found 
in Ref. 21. The explicit expressions for A a and B a can 
be cast into the form 



A„ 



B n = 



b a - ba + U(l - n ea - n oa ) 

t/2 (n (0) +n (0) )(2 _ n (0)_ n (0) )/4 



(14) 



(15) 



where n pa is the average number of impurity electrons 
in the corresponding channel p. It is calculated self- 
consistently by solving the equations 



'pa 



/(w)Im G pa {u). 
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n p °J are auxiliary numbers of particles defined as 



n<°> = - 



dw 



— /(u,)ImG<£>(u,) 



(17) 



with G p °J{uj) from Eq. $TTj). 
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; (18) 



4(n<?2 a + -1)EV 4 (4-^ 1 _ ff ) (0) 

(19) 

The superscript (0) in the correlation function in Eq. 
(fT9l) means that it is calculated at the Hartee-Fock level 
of approximation. The correlation functions in Eqs. ()18|) . 
(|19|) are obtained using the equation-of-motion method 
and the result is 

00 

E V k (a+_ CTCl _ CT (2ni_ CT - 1)) e* kRl = -Im f — 

— 00 

x m e _ A Gp _ H( 20) 

p— e,o ^ ' 

00 

E V k <a+^_ ff ) (0) e ikR > = -Im / — /(w) 

k J f 
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The coefficient A a [Eq. (fT4"|) ] ensures that the high- 
frequency limit of the exact Green's functions [Eq. (|7|)] 
coincides with the high-frequency limit of the corre- 
sponding approximate expressions [16[ . The approximate 
Green's functions have the correct limit of zero coupling 
between the impurities and the band electrons Vu — > 
(this limit is exactly solvable) [l6| . On the other hand, 
in the limit of strong coupling to the band electrons 
U/Vk — > 0, G pcr (uj) obviously coincides with the corre- 
sponding exact expression. Thus, the Green's functions 
we have obtained interpolate between the two exact lim- 
its of strong and weak coupling to the band electrons. 

In order to complete the construction of the impu- 
rity Green's functions we have to choose the parame- 
ter e. A discussion on the possible ways to fix its value 
in the case of the infinite-dimensional Hubbard model 
have shown that none of the possible choices is to be pre- 
ferred over the others [2l|. In a previous paper [r| on the 
two-impurity Anderson model (without spin-orbit inter- 
action), we have imposed the Friedel sum rule to deter- 
mine the value of e. The Friedel sum rule relates the aver- 
age impurity charge and the phase shift at the Fermi level 
at zero temperature and represents an exact condition on 
the to — > behavior of the Green's function. However, its 
application is confined only to the zero-temperature case 
while we intent to discuss the nonzero temperature be- 
havior of the model under consideration. Therefore, we 
choose the value of e from the condition n pa — n p °J . We 
have verified that, indeed, there is no significant differ- 
ence between the results obtained at T = with the two 
choices for e we have considered. Now the modified per- 
turbation theory for the TIAM with RSOI is completed. 
Let us stress that, unlike the studies of the two- impurity 
Kondo model, no degrees of freedom are integrated out 
and both the charge and the spin fluctuations are treated 
on equal footing in our approach. 

To proceed with the presentation of our results we need 
to specify the dispersion of the band electrons. The band 
electrons have a quadratic dispersion efe = k 2 /2m* — Eq 
where m* is the effective mass and Eq is the bottom of 
the band. The dispersion of the band electrons with given 
chirality can be cast into the form 

e kh = (k + hk R f/2m* -E -E R (22) 

with k R = m*a R being the Rashba wave vector and 
E R = k R /2m* is the Rashba energy. In the presence 
of the Rashba spin-orbit interaction the band bottom is 
actually at En + Eq. The Fermi wave vector is defined 
as fci? = \/2m*Ea (we choose the Fermi energy to be 0). 
We calculate the Hartree-Fock self-energy [Eq. ©] with 
the dispersion eth [Eq. (|22|) ]. The imaginary parts can 
be computed analytically but in order to obtain the real 
parts one has to resort to numerical integration. The re- 
sult for the imaginary part Im H P °J (w) is too complicated 
to be presented here. We shall only point out that it has 
a characteristic dependence on the band bottom of the 
form y/E R /{u + E + E R ) (see below) [12] . 



IV. RESULTS 

The numerical procedure involves the self-consistent 
solution of Eq. (|16[) . The impurity Green's functions de- 
pend on the values of the correlation functions in Eq. (f2TJ)) 
which in turn depend on G pa {ui). Therefore, another self- 
consistent procedure is required in order to compute the 
value of G p(J {uj) which is used in Eq. (fT6"j) . The procedure 
is rather effective and only in the case of small impurity 
separations kpR < 0.75 significant number of steps have 
to be performed in order to achieve self-consistency. In 
the following, we present our results for the set of param- 
eters e = -U/2, U = 3tt 7 , E q = 0.7C/ with 7 = V 2 m*/4. 
Note that the second-order perturbation theory in U is 
not expected to work for such large value of U. The 
choice for the value of e does not represent some spe- 
cial point for which a specific behavior can be expected 
(recall that the single-impurity Anderson model (SIAM) 
with e = —U/2 and with a symmetric density of states for 
the band electrons is a particle- hole symmetric model). 
We have pointed out previously [l6[ that for TIAM with 
e = —U/2 one obtains different physics compared to the 
case of symmetric SIAM. The reason is that in the case 
of TIAM there are two different effective energy levels in 
the even/odd channel given by e e / = e + Re S^ o (0) as 
well as two different elastic level widths 7 e / (see below). 
Thus, even with a symmetric density of states for the 
band electrons TIAM is not a particle-hole symmetric 
model. Moreover, in the present case, the 2D band elec- 
trons with the dispersion e^h [Eq. (|22|) ] are not actually 
a particle-hole symmetric system. These considerations 
suggest that the behavior of the TIAM with RSOI we 
are going to discuss will not qualitatively depend on any 
specific value of e. 

In Figs. 1 and 2, we show the spectral functions {i.e. 
— Im G p<7 (u))Iit) of the impurities' electrons in the even 
and odd channels at zero temperature T = 0. The results 
are presented in Fig. 1 for Rashba energy E R = 0.1J7 and 
for several impurity separations k F R and in Fig. 2 for 
impurity separation kpR = 1 and several values of E R . 
There are features in the left-hand side of both figures (at 
uj ~ —7.57 m Fig. 1) that result from the dependence 

of Im Y, P °J (uj) on the band bottom we have mentioned 
above. As evidenced from the figures, the height and the 
width of the Kondo peaks depend on both E R and k F R- 

The results for the dependence of the height and the 
width of the Kondo peaks on E R and kpR are presented 
in the insets of Fig. 2 and in Fig. 3. The k F R de- 
pendence shown in Fig. 3 is qualitatively similar to the 
dependence of the height and the width of the even/odd 
Kondo peaks on k F R for the 3D TIAM without RSOI 
[ill . The height of the even/odd Kondo peaks [Fig. 3, 
solid lines] increases/decreases with k F R for both values 
of E R and reaches extremum (maximum/minimum) for 
some value of k F R. Semi-quantitative understanding of 
this behavior has been suggested in Ref. 16. The relevant 
frequency range for the Kondo effect is the range of small 
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|u;| (\u>\ less than several times the Kondo temperature) 
and in it one can approximate the effective even/odd elas- 
tic level widths j p = — Im T,pJ (0) by 

{( 1+ i/ra)< 1± « 

+ ('-i/raJ' 1 **-'} ,23) 

where 

3+/ - = A ( M yf ± /T7f)). (24, 

In the case -Er = the effective level width is r y e / — 
27(1 ± Jo(kpR)). The extremums of this function are 
given by the zeros of J\{k F R) and the first zero is 
kpR = 3.8317 [22^ . Thus, 7 e has a minimum and 7 D has a 
maximum at this value of k F R. The data in Fig. 3 show 
that the height of the even/odd Kondo peaks reaches 
a maximum/minimum around k F R ~ 3.8 in agreement 
with the previous discussion When Er = 0.1J7 the 
minumum / maximum of 7 e / is realized for kpR 2.9 
and, indeed, the height of the even/odd Kondo peaks 
has a maximum/minimum in the vicinity of this value. 
Further, J {k F R) = for kpR = 2.4048 [22] and, con- 
sequently, for Er = 7 e = 7 for this value of fc^-R 
and both are equal to the elastic level width in the case 
kpR — > 00 which corresponds to the doubly-degenerate 
case of two identical infinitely-separated Anderson impu- 
rities. Fig. 3 shows that the height of the even and of 
the odd Kondo peak is almost equal to 1 in the vicinity 
of kpR s» 2.4 (recall that the data are normalized to the 
corresponding results for kpR 00 i.e. the case of a 
single Anderson impurity). 

The width of the even Kondo peak reaches maximum 
at k F R w 1.7 for En = and at somewhat smaller value 
of k F R sa 1.5 for Er = 0.1U [Fig. 3, lower panel, the 
dashed lines]. In a previous study [HI], we have related 
the occurrence of this maximum to the known change of 
the sign of the instantaneous spin-spin correlation func- 
tion which in the case of 3D TIAM without RSOI takes 
place at k F R — 37r/4 [13, EH- We are not aware of pre- 
vious works that have calculated the instantaneous spin- 
spin correlation function in the model we are studying. 
But, based on our previous experience, we can suppose 
that in the case of 2D TIAM with RSOI the instanta- 
neous spin-spin correlation function changes sign in the 
vicinity of k F R « 1.7 for En = and at k F R « 1.5 
for En = 0.1U being ferromagnetic/antiferromagnetic for 
lower /higher values of k F R. 

The width of the Kondo peak in the odd channel is a 
monotonically increasing function of k F R [Fig. 3, upper- 
panel, the dashed lines] . It has a maximum at the same 
value of k F R at which a minimum of the height of the odd 
Kondo peak is realized (see above) . At the same value of 
k F R there is a maximum of the height and a minimum 



of the width of the even Kondo peak [Fig. 3]. Physi- 
cally, maximum of the height and minimum of the width 
of the even Kondo peak means that it is most difficult 
for the even Kondo resonance to be destroyed by the an- 
tifcrromagnetic correlations of the impurity spins at the 
corresponding impurity separation. On the other hand, 
the odd Kondo peak has minimal height and maximal 
width for the same impurity separation, that is, the an- 
tiferromagnetic impurity correlations tend to most easily 
destroy the odd Kondo resonance. 

The results for the dependence of the height and the 
width of the even/odd Kondo peaks on the Rashba en- 
ergy En are presented in the insets of Fig. 2. The En de- 
pendence is qualitatively the same in both the even and 
the odd channel. The height has a maximum at some 
value of En (different for the even and the odd channel) 
and the width is a monotonically increasing function of 
En- 

Finally, we show the results for the dependence of the 
Kondo temperature Tk on En [Fig. 4]. The Kondo 
temperature has been determined from the variation of 
the width of the Kondo peak in the even/odd channel 
with the temperature. We are not aware of previous 
works that have considered the reliability of obtaining 
Tk in the framework of the modified perturbation the- 
ory. Nevertheless, we think that this approach provides 
reasonable predictions for the dependence of Tk on the 
model parameters. As evidenced from Fig. 4, Tk has a 
non-monotonic dependence on En for sufficiently small 
impurity separation k F R < 1.2 — 1.3. For so small im- 
purity separation there is a significant difference between 
the Kondo temperatures in the even and in the odd chan- 
nels with Tkc > Tko- Recall that in the limit k F R — > 
(the two impurities merge at some point) the odd chan- 
nel disappears and the even channel corresponds to one 
impurity with elastic level width j e (oj). For larger val- 
ues of Er > 0.1 — 0.12U the dependence of Tk is almost 
linear. For larger impurity separations the even and the 
odd Kondo temperature are practically equal and their 
En dependence is monotonic - Tk increases almost lin- 
early with En- Thus, we predict that the inclusion of 
the RSOI will not lead to an exponential increase of Tk- 
For larger values of En the Kondo temperature increases 
almost linearly for any value of R similarly to the result 
obtained in a different approach in the single-impurity 
Anderson model [l2j]. 

To conclude, we have constructed a modified pertur- 
bation theory for the 2D two-impurity Anderson model 
with Rashba spin-orbit interaction. The ansatz interact- 
ing self-energy is chosen in such a way that the impurity 
Green's functions are exact up to order w -4 . The height 
and width of the Kondo peaks in the even/odd chan- 
nels are obtained numerically as functions of the inter- 
impurity distance and the strength of the Rashba spin- 
orbit interaction. It is predicted that the Kondo tem- 
perature will not have an exponential increase with the 
Rashba energy instead having only almost linear depen- 
dence on it. 
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FIG. 1: (Color online) The impurity spectral functions in the 
even/odd channels for Er — 0.1(7 and several values of UfR. 
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FIG. 2: (Color online) The impurity spectral functions in the 
even/odd channels for UfR = 1 and several values of Er. The 
insets: The Eh dependence of the height (solid line) and the 
width (dashed line) of the even/odd Kondo peaks. 
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FIG. 3: (Color online) The UfR dependence of the height 
(solid line) and the width (dashed line) of the even/odd Kondo 
peaks for Er = and Er — 0.1)7. The data are normal- 
ized with respect to the corresponding quantities in the limit 
kpR —¥ oo. 
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FIG. 4: The Er dependence of the Kondo temperature for 
two values IcfR. The data are normalized with respect to the 
Kondo temperature for Er — in the limit fc^-R — > oo. 



